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SUMMARY 

The  basis  for  this  research  is  that  rockbursts  produce  seismic  signals  that  are 
similar  to  those  associated  with  near  surface  earthquakes  and  explosions.  A  complete 
understanding  of  source  mechanisms  is  necessary  for  the  correct  identification  of  these 
events.  Some  of  the  methods  developed  for  source  mechanism  studies  also  encompass 
imaging  techniques  appropriate  for  nuclear  event  verification,  and,  as  a  result,  have 
implications  for  the  Comprehensive  Test  Ban  Treaty  (CTBT). 

PCI/Atlantic  Scientific  INSAR  software  was  purchased  last  fall  for  elevation 
change  detection  at  a  rockburst  site  using  radar  interferometry.  We  have  been 
familiarizing  ourselves  with  this  software  and  radar  remote  sensing.  Images  have 
recently  been  acquired  that  cover  the  site  of  a  large  rockburst. 

The  air  wave  can  interfere  with  the  rock  wave  in  Ground  Penetrating  Radar  (GPR) 
transillumination  experiments  making  it  difficult  to  pick  the  rock  wave  first  arrivals  in 
mines.  This  was  particularly  a  problem  at  small  transmitter  to  receiver  distances.  The 
distinction  between  these  arrivals  could  be  overcome  by  the  use  of  higher  frequency 
antennas  at  small  transmitter  to  receiver  distances.  Lower  frequency  antennas  might  still 
be  necessary  for  longer  transmitter  to  receiver  distances  to  obtain  the  necessary  depth  of 
penetration.  A  comparison  of  the  tomograms  from  the  1997  and  1998  GPR  experiments 
at  the  Stillwater  Mine,  Montana,  will  help  us  to  decide  if  we  can  use  this  method  to  detect 
change  in  a  pillar.  A  comparison  with  electrical  measurements  on  the  same  pillar  will 
help  with  this  interpretation. 

Past  experiments  have  demonstrated  a  clear  relationship  between  Barkhausen 
Noise  (BN)  and  stress  in  pipeline  steel  samples.  Stress  can  affect  the  nature  of  the 
pinning  sites  at  domain  boundaries  which  in  turn  plays  a  role  in  the  way  that  domain 
walls  move  and  cause  BN  pulses.  We  have  conducted  experiments  to  detect  a 
relationship  between  BN  and  stress  in  magnetite-rich  rocks.  These  experiments  were 
performed  by  applying  uniaxial  loading  from  4500N  to  55000N  to  cylindrical  core 
samples  of  magnetite-rich  rocks  obtained  from  the  Stillwater  Complex  in  Montana.  A 
pick-up  coil  was  used  to  detect  the  BN  pulses  and  an  electromagnet  was  used  to 
magnetize  the  samples.  BN  pulses  were  recorded  during  a  specified  time  interval  for 
each  increment  of  loading.  The  pick-up  coil  was  positioned  to  be  near  magnetite-rich 
regions  of  the  sample.  We  have  observed  a  linear  increase  in  the  number  of  BN  pulses 
with  increased  loading.  Ongoing  research  is  focused  on  refining  the  experimental 
procedure  to  increase  the  bandwidth  of  recording  and  eliminate  possible  sources  of  noise 
in  the  experiments.  The  results  of  this  study  might  help  us  to  understand  the  nature  of  the 
remagnetization  that  often  accompanies  large-scale  tectonic  events.  Moreover,  this  study 
could  lead  to  a  method  of  nondestructive  testing  to  determine  state  of  stress  at  mine  walls. 

A  series  of  electrical  measurements  have  been  made  on  a  new  pillar  of  rock  in  the 
Stillwater  Mine  over  a  period  of  two  years.  The  measurements  have  been  used  to 
generate  models  of  the  resistivity  in  the  pillar  and  the  surrounding  rock. 

The  models  show  a  cyclic  or  seasonal  change  in  resistivity  of  some  of  the  rock 
surrounding  the  pillar.  This  was  unexpected  ,  given  the  depth  of  the  rock  below  the 
surface.  The  resistivity  in  the  interior  of  the  pillar  showed  a  general  decrease  over  time. 
The  decrease  is  attributed  to  the  opening  of  cracks  and  the  redistribution  of  pore  water. 
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Soiu-ce  mechanisms  for  an  addition  15  rockbursts  in  the  Lucky  Friday  Mine  of 
Northern  Idaho  have  been  determined.  The  rockburst  events  had  equivalent  Richter 
magnitudes  ranging  from  approximately  2.5  to  3.5,  and  in  some  cases  caused  extensive 
damage  to  mine  workings.  All  of  the  events  were  within  5700  ft  of  the  surface,  and  would 
therefore  be  classified  as  extremely  near  surface  seismic  events.  Naturally  occurring 
earthquakes  in  this  geographic  region  have  focal  mechanisms  much  deeper  than  the  rockburst 
events.  All  of  the  rockburst  source  mechanisms  have  equivalent  double  couple  sources  similar 
to  those  caused  by  earthquakes.  None  of  the  rockburst  events  can  be  explained  by  implosive 
or  explosive  types  of  source  mechanisms. 

Identification  of  the  state  of  rock  bolts  inside  boreholes  is  a  difficult  problem. 

Rock  bolts  that  are  not  seated  properly  inside  a  borehole  can  fail  imder  a  variety  of 
conditions  resulting  in  catastrophic  cave-ins  and  rockbursts.  For  this  study,  modeled 
frequency  response  curves  were  compared  to  closed-form  solutions  for  axial  vibrations, 
lateral  vibrations,  and  strings  clamped  at  both  ends.  In  addition,  the  computer  results 
were  compared  to  results  from  a  laboratory  experiment  utilizing  two  different  rock  bolts. 
Finally,  rock  bolts  from  the  Stillwater  Mine  were  modeled  and  results  were  obtained  for 
two  seating  situations. 

Using  a  Tinius-Ohlsen  testing  machine,  a  Simmons  eye  rock  bolt  and  DYWIDAG 
rock  bolt  were  placed  under  varying  amounts  of  tension  to  investigate  frequency  mode 
response.  Impulse  forces  were  applied  as  hammer  taps  to  the  ends  and  sides  of  the  bolts  and 
the  resulting  sounds  were  recorded  with  an  audiocassette  recorder.  These  recordings  were 
analyzed  in  the  frequency  domain  and  compared  with  computer  modeling  results. 

Investigation  focused  on  axial  responses  rather  than  bending  or  transverse  responses 
due  to  higher  signal-to-noise  ratio  for  axial  impulses.  Frequency  peaks  occurred  at  4200  Hz 
and  3450  Hz  for  the  Simmons  eye  bolt  axial  modes.  The  effect  of  clamp  placement  on  the 
bolt  and  non-direct  hits  on  the  bolt  end  most  likely  contribute  to  the  combination  of  axial  and 
bending  peaks. 

Impulse  response  measurements  for  the  DYWIDAG  rock  bolt  showed  frequency 
peaks  at  750  Hz,  1250  Hz,  and  2100  Hz.  Modeling  indicated  that  the  combination  of  rock- 
bolt  ridges  and  clamp  placement  contributed  to  the  multiple  peaks  and  lowered  the 
frequencies. 

In  addition  to  laboratory  measurements,  hammer  impulse  tests  were  performed  on  in 
situ  rock  bolts  in  the  Stillwater  mine.  Underground  miners  maintain  that  they  can  predict  mine 
wall  loading  by  the  sound  of  a  rock  hammer  hitting  a  rock  bolt.  The  fundamental  frequency 
contained  in  the  mine  recordings  was  1260  Hz.  This  frequency  peak  occurred  within  a  400 
Hz  range.  Different  rock-bolt  lengths,  ridges  on  the  sides  of  the  rock  bolts,  non-direct  hammer 
strikes  on  the  ends  of  the  rock  bolts,  and  variable  coupling  of  the  bolts  to  the  rock  all 
contribute  to  the  variation  in  frequencies.  Due  to  the  variable  nature  of  rock-bolt  coupling,  the 
range  in  frequencies  observed  from  the  mine  rock  bolts  is  probably  not  attributable  to  changes 
in  loading. 
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REMOTE  SENSING  METHODS  APPLIED  TO  ROCKBURST  MONITORING 


Bacl^round 

Remote  sensing  methods  were  to  be  employed  in  order  to  observe  any  changes 
that  might  occur  at  the  surface  above  a  mining  area  as  a  result  of  a  rockburst.  Originally, 
hyperspectral  or  multispectral  images  were  to  be  purchased  for  temporal  change  detection 
over  a  rockburst  site;  however,  due  to  an  early  budget  reduction  in  the  proposal,  there  was 
not  enough  money  allocated  to  carry  out  this  part  of  the  project.  In  addition,  it  became 
evident  early  on  that  rather  than  focusing  on  spectral  change  detection,  it  might  be  more 
advantageous  to  look  at  surface  displacements  using  the  relatively  new  technique  of  radar 
interferometry.  Consequently,  the  time  allocated  to  this  project  was  spent  familiarizing 
us  with  radar  interferometry  and  trying  to  secure  funds  from  another  source  in  order  to 
carry  out  the  remote  sensing  part  of  the  project. 

Update 

At  the  beginning  of  this  report  period,  NASA  funds  were  obtained  through  a 
JOVE  Augmentation  Grant  to  help  support  this  part  of  this  project.  As  a  consequence, 
PCI/Atlantic  Scientific  INSAR  software  was  purchased  for  elevation  change  detection  at 
a  rockburst  site  using  radar  interferometry.  We  have  been  familiarizing  ourselves  with 
this  software  and  radar  remote  sensing.  Images  have  recently  been  acquired  that  cover 
the  site  of  a  large  rockburst. 
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GROUND  PENETRATING  RADAR  MONITORING 


Introduction 

We  conducted  Ground  Penetrating  Radar  (GPR)  tomography  experiments  at  the 
Stillwater  Mine  in  Montana  to  see  if  we  could  detect  pillar  changes  using  this  method. 
Experiments  were  completed  on  July  8  and  9,  1997,  January  13,  1998,  and  August  26, 
1998. 

The  Stillwater  mine  is  a  platinum-palladium  underground  mine  located  at  the 
edge  of  the  Beartooth  Mountains,  near  Nye,  Montana.  This  mine  is  the  only  platinum- 
palladium  mine  in  operation  in  the  United  States.  Here,  platinum  and  palladium  are 
mined  from  the  Stillwater  Complex.  The  rocks  of  this  large,  layered,  mafic  intrusive 
complex  are  predominantly  gabbro  and  related  rock  types.  This  igneous  intrusion  runs 
along  the  northwest  margin  of  the  Beartooth  Mountains  in  Stillwater  and  Sweetgrass 
counties.  The  middle  Cambrian  intrusion  exhibits  a  west  to  northwest  strike  with  a  50  to 
90  degree  dip.  The  complex  is  host  to  pyroxenes,  gabbros,  norites,  diorites  and 
anorthosites  intruded  into  Precambrian  rocks  and  overlain  by  Paleozoic  sediments  (Yates, 

1984) .  Historically,  this  mine  has  not  been  prone  to  rockbursts;  but,  as  mining  has 
•progressed  deeper,  rock  failure  has  become  more  of  a  threat. 

Procedure 

In  these  tests,  the  transillumination  mode  of  data  gathering  was  employed  in 
which  the  transmitting  and  receiving  antennas  were  placed  on  opposite  sides  of  a  pillar 
(Figure  1).  The  receiving  antenna  remained  fixed  while  the  transmitter  was  moved  to 
create  a  total  of  fifteen  receiver  gathers.  This  method  of  data  gathering  proved  to  be 
efficient.  Table  1  shows  the  field  recording  parameters  used  for  this  study. 

Arrival  times  of  the  rock  wave  were  picked  by  hand  from  the  receiver  gathers. 
These  times,  together  with  the  source  and  receiver  locations,  will  be  used  to  invert  for 
velocity  by  the  Simultaneous  Iterative  Reconstruction  Technique  (SIRT)  (Peterson  et  al., 

1 985) .  Software  developed  at  the  former  United  States  Bureau  of  Mines  was  used  to 
invert  for  velocity  structure  of  the  pillar  for  the  1997  experiment.  The  software  is  called 
3DTOM  and  runs  under  DOS  on  personnel  computers  (Jackson  and  Tweeton,  1996).  A 
graduate  student  is  processing  and  interpreting  these  data  as  part  of  a  thesis  projeet. 
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Table  1 .  Data  Collection  Parameters 


Samples 

1024 

Sample  Frequency 

1008.94  MHz 

Stacks 

128 

Antennas 

100  MHz 

Antenna  orientation 

Vertical 

The  velocity  tomogram  for  the  1997  experiment  is  shown  in  Figure  2  and  the 
corresponding  raypath  density  plot  is  shown  in  Figure  3.  These  data  were  processed 
using  SIRT  with  a  straight  raypath  forward  model.  The  termination  criteria  was  15 
iterations  or  a  residual  tolerance  of  0.01,  whichever  occurs  first.  A  nodal  constraint  of 
300  m/ps  was  used  for  nodes  occurring  in  the  tunnel  region  adjacent  to  the  pillar  face. 
This  is  the  electromagnetic  wave  speed  of  free  space. 

The  tomogram  (Figure  2)  shows  the  outline  of  the  pillar  with  higher  velocities 
occurring  near  the  pillar  face.  Ordinarily,  we  would  observe  high  relative  velocities  at  the 
nose  of  the  pillar  due  to  increased  pore  volume  caused  by  fracturing.  Low  velocities, 
however,  are  observed  at  the  nose  of  the  pillar  in  the  lower  left  of  Figure  2.  This  is  a 
correspondingly  area  of  low  raypath  coverage  (Figure  3)  and  this  might  be  an  artifact  of 
the  inversion. 
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Fig.  2.  GPR  wave-speed  tomogram  for  pillar  study  at  Stillwater  Mine. 
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Fig.  3.  Raypath  sample  density  plot.  Contours  show  the  outline  of  the 
pillar  in  the  middle  of  a  rectangular  model  grid. 


Discussion  and  Conclusions 

Reflections  from  the  tunnel  walls  as  well  as  scattering  from  metal  in  the  tunnel 
can  often  interfere  with  reflections  received  fi-om  discontinuities  in  the  rock.  Preliminary 
tests  of  a  wire  mesh  (“chicken  wire”)  to  help  shield  our  antennas  proved  inconclusive. 

The  air  wave  can  interfere  with  the  rock  wave  in  transillumination  experiments 
making  it  difficult  to  pick  the  rock  wave  first  arrivals.  This  was  particularly  a  problem  at 
small  transmitter  to  receiver  distances.  The  distinction  between  these  arrivals  could  be 
partly  overcome  by  the  use  of  higher  frequency  antennas  at  small  transmitter  to  receiver 
distances.  Lower  frequency  antennas  might  still  be  necessary  for  longer  transmitter  to 
receiver  distances  to  obtain  the  necessary  depth  of  penetration.  Due  to  the  limited  depth 
of  penetration  observed  in  the  Stillwater  Mine  we  must  use  100  MHz  antennas  for  this 
pillar. 

A  comparison  of  the  tomograms  from  the  1997  and  1998  GPR  experiments  will 
help  us  to  decide  if  we  can  use  this  method  to  detect  changes  in  the  pillar.  A  comparison 
with  electrical  measurements  on  the  same  pillar  will  help  with  this  interpretation. 
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BARKHAUSEN  NOISE  MONITORING 
Introduction 

Barkhausen  Noise  (BN)  refers  to  the  electromagnetic  pulses  produced  when 
magnetic  domain  walls  move  in  a  magnetic  material  in  response  to  the  applied  magnetic 
field.  Professor  Heinrieh  Barkhausen  (1919)  demonstrated  that  the  magnetization  of  a 
sample,  which  is  characterized  by  a  hysteresis  curve,  is  not  continuous  but  in  fact  is  made 
up  of  a  series  of  steps  which  are  produced  as  domain  walls  move  in  response  to  the 
applied  magnetic  field  (Figure  4).  We  have  conducted  experiments  to  detect  BN  in 
magnetite-rich  rock  subjected  to  loading. 
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Fig.  4.  Generalized  hysteresis  eurve  for  a  magnetie  material  showing 
Barkhausen  noise  (B — magnetic  induction;  H — magnetic  field  strength; 
after  Bozorth,  1951). 


Domains 

Domains  have  their  own  magnetic  strength  and  direction.  The  sum  of  the  strengths 
and  direetions  of  all  the  domains  results  in  the  total  magnetization  of  the  entire  sample. 
Generally,  the  domains  are  different  sizes  and  have  different  directions  of  magnetization, 
which  is  typical  for  ferromagnetic  and  ferrimagnetic  material.  Domain  walls,  shown  in  black, 
in  Figure  5,  separate  the  individual  domains  from  each  other.  Domains  walls  become  fixed  or 
piimed  at  crystalline  imperfeetions.  When  a  ferromagnetic  material  is  placed  in  a  magnetic 
field,  the  domain  walls  move.  The  movement  of  the  domain  walls  increases  the  size  of  the 
domains  that  have  the  same  magnetization  direction  as  the  magnetic  field,  and  decreases  the 
size  of  the  domains  with  opposing  directions  (Figure  5b). 
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Fig  5.  (a)  Example  of  unmagnetized  domains.  The  direction  of  the  arrows 
indicates  the  direction  of  the  magnetization  of  the  sample.  The  black  areas 
are  the  domain  walls,  which  move  when  the  material  is  magnetized,  (b) 

Domains  magnetized  by  the  field  H.  The  fine  arrows  indicate  the  direction 
of  the  motion  of  the  domain  walls. 

This  movement  increases  the  overall  magnetization  of  the  field  in  the  direction  of 
H,  giving  rise  to  a  hysteresis  curve.  The  domain  wall  movement,  however,  is  not  smooth. 
A  wall  moves  until  it  reaches  a  pinning  site,  at  which  point  it  is  held  up  until  the  magnetic 
field  is  large  enough  to  move  the  wall  past  the  site.  Once  the  site  is  overcome  the  wall 
moves  smoothly  imtil  it  encounters  another  pinning  site.  This  process  continues  until  the 
sample  is  saturated.  These  quick  wall  movements,  or  jumps  over  the  pinning  sites,  in 
ferromagnetic  material  can  induce  an  electromagnetic  pulse  called  BN. 

Barkhausen  Noise  and  Stress 

BN  has  been  related  to  grain  size  (Tiitto,  1977;  Ranjan  et  al.,  1987),  carbon 
content  (Ranjan  et  al.,  1987),  pearlite  colony  and  lamellar  spacing  (Svmdstrom  and 
Torronen,  1979;  Clapham  et  al.,  1991),  and  residual  stress  (Sandstrom  and  Torronen, 
1979;  Dhar  and  Atherton,  1993).  Elastic  stresses  influence  how  domains  choose  their 
easy  direction  of  magnetization.  As  the  result  of  this  magnetoelastic  interaction, 
compressive  stresses  will  decrease  the  intensity  of  BN  and  tensile  stresses  increase  the 
intensity  of  BN  (Jagadish  et  al.,  1990).  Consequently,  a  measurement  of  the  intensity  of 
BN  in  a  sample  can  provide  information  about  the  residual  stress.  BN  is  commonly 
observed  by  using  a  pick-up  coil  and  an  amplifier  (Figure  6). 
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Fig.  6.  Method  of  observing  Barkhausen  noise  (after  Bozorth,  1951). 


Magnetoelastic  testing  using  BN  measurements  have  been  used  primarily  in  the 
nondestructive  testing  of  pipeline  steel  (Tiitto,  1991;  Jagadish  et  ah,  1990;  Krause  et  al., 
1995;  Mandel  et  ah,  1996).  This  method  of  nondestructive  testing  has  not  been  tried  on 
earth  materials,  although  the  Barkhausen  effect  has  been  observed  in  rocks.  Furthermore, 
some  have  suggested  that  domain  wall  adjustments  bring  about  remagnetization  during 
tectonic  events  in  response  to  the  changing  stress  field  (Hudson  et  ah,  1989;  Borradaile, 
1996). 

Experimental  Procedure 

A  block  diagram  of  the  circuit  that  was  built  to  record  BN  can  be  seen  in  Figure  7. 
An  electromagnet  with  a  10  Hz  driving  signal  was  used  to  generate  a  magnetic  field  of 
about  86,000  gammas.  A  read  head  from  a  Sony  Walkman  was  used  to  detect  BN.  The 
data  were  recorded  using  a  computer  sound  card  with  a  sample  rate  of  22  KHz.  Sample 
data  and  the  10  Hz  driving  signal  are  shown  in  Figure  8.  Figure  9  shows  filtered  data 
with  BN. 

A  compressional  load  was  applied  to  a  cylindrical  core  samples  that  were  obtained 
from  the  Stillwater  Complex  in  Montana.  The  core  samples  were  selected  to  be 
magnetite  rich.  Figure  10  shows  the  experimental  setup.  A  Tinius-Olsen  Hydraulic 
Tension  and  Compression  Tester  was  used  to  load  the  sample.  The  read  head  was 
oriented  to  be  near  magnetite-rich  portions  of  the  samples.  BN  was  recorded  for  a  small 
time  interval  at  each  1000  lb.  increment  of  loading. 

Results  and  Conclusions 

Figure  1 1  shows  the  total  number  of  pulses  recorded  at  each  load  level  as  a 
function  of  load.  We  have  observed  a  linear  increase  in  the  number  of  BN  pulses  with 
increased  loading.  Ongoing  research  is  focused  on  refining  the  experimental  procedure  to 
increase  the  bandwidth  of  recording  and  eliminate  possible  sources  of  noise  in  the 
experiments.  The  results  of  this  study  might  help  us  to  understand  the  nature  of  the 
remagnetization  that  often  accompanies  large-scale  tectonic  events.  Moreover,  this  study 
could  lead  to  a  method  of  nondestructive  testing  to  determine  state  of  stress  at  mine  walls. 
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Frequency  Generator  Amplifier  Electromagnet 


Fig.  7.  Experimental  setup  to  detect  BN. 


Fig.  8.  Unfiltered  record  showing  the  higher  amplitude  driving  signal  and 
the  signal  recorded  from  the  read  head. 


Number  of  Samples  =  1  second 


Number  of  Samples  =  1  second 


Number  of  Samples  =  1  second 


Fig.  9.  Three  tests  with  BN  pulses,  (a)  Magnetized  read  head  only  (b)  A  nail 
placed  in  the  magnetic  field,  (c)  A  core  sample  containing  magnetite  in  the 
magnetic  field. 
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Uniaxial  Loader 


Fig.  10.  Experimental  loading  setup. 


Total  Number  of  Pulses 


1350 


Fig.  1 1 .  Total  number  of  BN  pulses  for  each  load  level  (•).  A  least  squares  best  fit  line 
has  been  added  to  show  an  increasing  linear  trend  in  the  data. 
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ELECTRICAL  RESISTIVITY  MONITORING 
Introduction 

The  electrical  resistivity  monitoring  of  a  pillar  of  rock  on  the  3200  level  of  the 
Stillwater  Mine,  Nye,  MT  has  been  in  progress  for  two  years.  During  that  time,  seven  sets 
of  measurements  have  been  made.  The  pillar  is  at  depth  of  about  2000  feet  below  the 
ground  surface  and  was  excavated  shortly  before  the  observations  began.  The  pillar  is 
roughly  triangular  in  cross  section.  (See  Figure  1.)  Eight  electrodes  were  utilized  to  make 
electrical  measurements  of  the  pillar  and  surrounding  rock.  These  measurements  were 
used  to  develop  a  tomographic  model  of  the  pillar  and  the  neighboring  rock.  The  model  is 
a  region  15  by  18  meters,  subdivided  into  270  blocks,  1  meter  on  a  side. 

Results 

Figure  12  shows  a  surface  plot  of  the  resistivity  (vertical  axis)  versus  the 
horizontal  position  in  the  model.  The  electrode  positions  are  shown  as  circles  and  the 
outline  of  the  pillar  is  shown  by  the  line  connecting  the  circles.  The  model  shows  high 
resistivity  (3,000  ohm-m)  in  the  interior  of  the  pillar,  low  resistivity  (200  ohm-m)  in  the 
timnels  alongside  the  pillar  and  intermediate  resistivity  (1,000  ohm-m)  in  the  surrounding 
rock. 

Over  the  first  six  months  (Figures  12  to  1 5)  there  is  a  gradual  increase  in  the 
resistivity  of  the  rock  outside  the  pillar  from  1,000  ohm-m  to  3,000  ohm-m  and  some 
increases  and  decreases  within  the  pillar.  By  the  time  one  year  has  passed,  (Figure  16)  the 
region  exterior  to  the  pillar  has  returned  to  values  similar  to  that  at  the  start  (Figure  12) 
and  the  values  within  the  pillar  are  somewhat  lower. 

During  the  second  year  (Figures  16  to  18),  the  resistivity  in  the  region  outside  the 
pillar  increases  in  the  winter  (Figure  17)  and  decreases  in  the  summer  (Figure  1 8),  similar 
to  the  cycle  in  the  first  year.  This  cyclic  behavior  in  the  exterior  region  suggests  a  source 
connected  with  the  weather  conditions  at  the  surface.  One  possibility  is  that  the  freezing 
of  the  precipitation  and  the  shallow  groundwater  removes  a  source  of  porewater  at  depth 
(2,000  feet),  resulting  in  a  increase  in  the  resistivity.  The  apparent  in-phase  response  of 
the  rocks  at  depth  to  the  available  liquid  water  at  the  surface  implies  a  well  connected 
hydraulic  system. 

The  resistivity  changes  in  the  pillar  are  more  easily  seen  by  plotting  the  ratio  of 
the  modeled  resistivity  at  a  particular  date  to  the  resistivity  at  the  start  of  the  series.  These 
results  are  shown  in  Figures  19  to  24  as  contour  plots  of  log2(resistivity  ratio).  An 
increase  (decrease)  of  one  unit  is  then  equal  to  an  increase  of  2x  (decrease  of  ‘Ax  )  in  the 
resistivity.  Figures  19, 20  and  21  show  a  pattern  of  decreasing  resistivity  in  the  central 
part  of  the  pillar,  while  at  the  same  time  the  resistivity  is  increasing  outside  of  the  pillar. 
The  decrease  is  still  present  in  the  summer  of  1998  (Figure  22)  and  persists  the  following 
winter  (Figure  23).  The  effect  is  present  the  next  summer  with  some  decrease  in 
magnitude  (Figure  24)  and  there  is  the  suggestion  that  the  low  is  moving  from  the  pillar 
into  the  surrounding  rock. 
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Discussion  and  Conclusions 

The  monitoring  and  modeling  of  pillar  region  for  two  years  has  revealed  an 
apparent  seasonal  cycle  in  some  of  the  rock  surrounding  the  pillar  and  a  general  decrease 
in  the  resistivity  of  the  pillar.  There  is  a  possibility  that  this  decrease  is  migrating  from 
the  pillar  to  the  surrounding  rock. 

In  the  last  report,  a  discussion  of  some  other  observed  resistivity  changes  in 
response  to  stress  did  not  seem  to  help  in  the  interpretation  of  this  data  from  the  Stillwater 
Mine.  Temperature  effects  this  deep  in  a  mine  also  seem  unlikely.  It  was  suggested  that 
opening  of  cracks  and  redistribution  of  pore  water  might  be  a  more  plausible  explanation. 
The  apparent  response  of  some  of  the  surrounding  rock  to  the  availability  of  surface  water 
was  completely  unexpected. 
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Fig.  12.  Surface  plot  of  resistivity  in  a  pillar  in  the  Stillwater  Mine, 
positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown 
the  circles. 
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Fig.  13.  Surface  plot  of  resistivity  in  a  pillar  in  the  Stillwater  Mine.  The  electrode 
positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line  connecting 
the  circles. 
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Fig.  15.  Surface  plot  of  resistivity  in  a  pillar  in  the  Stillwater  Mine.  The  electrode 
positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line  connecting 
the  circles. 
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Fig.  16.  Surface  plot  of  resistivity  in  a  pillar  in  the  Stillwater  Mine.  The  electrode 
positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line  connecting 
the  circles. 
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Fig.  17.  Surface  plot  of  resistivity  in  a  pillar  in  the  Stillwater  Mine.  The  electrode 
positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line  connecting 
the  circles. 
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Fig.  19.  Contour  plot  of  log2(resistivity  ratio)  for  a  pillar  in  the  Stillwater  Mine.  The 
electrode  positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line 
connecting  the  circles. 
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Fig.  20.  Contour  plot  of  log2(resistivity  ratio)  for  a  pillar  in  the  Stillwater  Mine.  The 
electrode  positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line 
connecting  the  circles. 
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Fig.  21.  Contour  plot  of  log2(resistivity  ratio)  for  a  pillar  in  the  Stillwater  Mine.  The 
electrode  positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line 
connecting  the  circles. 
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Fig.  22.  Contour  plot  of  log2(resistivity  ratio)  for  a  pillar  in  the  Stillwater  Mine.  The 
electrode  positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line 
connecting  the  circles. 
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Fig.  23.  Contour  plot  of  log2(resistivity  ratio)  for  a  pillar  in  the  Stillwater  Mine.  The 
electrode  positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line 
connecting  the  circles. 
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Fig.  24.  Contour  plot  of  log2(resistivity  ratio)  for  a  pillar  in  the  Stillwater  Mine.  The 
electrode  positions  are  shown  as  circles  and  the  outline  of  the  pillar  is  shown  by  the  line 
connecting  the  circles. 
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ROCKBURST  SOURCE  MECHANISM  STUDIES 
Summary 

Source  mechanisms  for  an  addition  15  rockbursts  in  the  Lucky  Friday  Mine  of 
Northern  Idaho  have  been  determined.  The  rockburst  events  had  equivalent  Richter 
magnitudes  ranging  from  approximately  2.5  to  3.5,  and  in  some  cases  caused  extensive 
damage  to  mine  workings.  All  of  the  events  were  within  5700  ft  of  the  surface,  and  would 
therefore  be  classified  as  extremely  near  surface  seismic  events.  Naturally  occurring 
earthquakes  in  this  geographic  region  have  focal  mechanisms  much  deeper  than  the  rockburst 
events.  All  of  the  rockburst  source  mechanisms  have  equivalent  double  couple  sources  similar 
to  those  caused  by  earthquakes.  None  of  the  rockburst  events  can  be  explained  by  implosive 
or  explosive  types  of  source  mechanisms. 

Background 

An  array  of  geophones  is  distributed  throughout  the  Lucky  Friday  Mine  near  Mullen, 
Idaho.  The  positions  of  the  geophones  are  restricted  by  the  accessibility  of  mine  openings,  so 
the  spatial  distribution  of  the  detectors  is  limited.  Researchers  at  the  U.  S.  Bureau  of  Mines 
have  recognized  the  problems  associated  with  using  geophones  within  mines  for  source 
mechanism  studies.  For  source  mechanism  determinations,  not  all  quadrants  are  equally 
covered,  and  optimal  distances  may  not  be  obtained.  Some  of  the  geophones  are  near  the 
elevation  of  the  rockbursts,  and  the  result  is  that  there  are  uncertainties  in  stereonet 
projections.  There  are  often  several  geophones  near  nodal  planes  for  various  seismic  wave 
phases,  so  that  small  location  errors  in  the  source  position  cause  large  changes  in  source 
mechanism  estimation.  In  the  case  of  the  Lucky  Friday  Mine,  there  is  an  independent  seismic 
system  that  has  many  geophones  located  near  the  working  faces  of  the  mine,  £ind  that  system 
is  utilized  to  locate  events  within  a  few  meters.  The  independent  seismic  system  is  called  the 
micro-seismic  system,  and  the  source  mechanism  system  is  called  the  macro-seismic  system. 
In  most  cases,  the  effects  of  source  location  error  may  be  minimized. 

Results 

Table  2  lists  studied  rockburst  events.  All  of  the  rockburst  events  occurred  within 
5500  ft  of  the  surface.  Therefore  all  of  the  mining  induced  seismic  events  tend  to  have 
shallower  epicenters  than  the  natural  earthquakes  of  the  region.  All  of  the  rockbursts  have 
source  mechanisms  that  are  similar  to  standard  double  couple  events  produced  by  tectonic 
activity,  and  therefore  the  source  mechanisms  are  similar  to  earthquake  sources.  None  of  the 
events  appear  to  have  either  implosive  or  explosive  source  mechanisms.  Standard  “beach  ball’ 
type  illustrations  used  to  show  earthquake  events  are  shown  in  Appendix  A,  and  are  used  to 
show  the  possible  double  couple  mechanisms  for  each  rockburst  event  studied.  In  some  cases, 
more  than  one  solution  is  possible.  When  multiple  solutions  are  found,  they  are  included.  The 
first  solution  shown  for  each  event  is  the  preferred  solution  because  that  particular  one  is  most 
compatible  with  the  geometry  of  the  mine  workings. 
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Table  2.  Table  of  Events 


DATE 

TIME 

MAGNITUDE 

COMMENTS 

1-13-95 

12:36:07 

Dip  slip 

AM 

N/S 

5-27-95 

04:17:22 

Dip  slip 

AM 

■ 

E/W 

6-22-95 

09:38:18 

Dip  slip 

PM 

E/W 

9-19-95 

01:40:37 

Strike  slip 

PM 

NW/SE 

2-8-96 

06:01:30 

3.5 

Strike  slip 

PM 

NE/SW 

3-2-96 

12:37:35 

2.6 

Dip  slip 

AM 

NE/SW 

4-2-96 

10:00:01 

2.8 

Strike  slip 

PM 

NW/SE 

5-30-96 

03:00:40 

3.5 

Strike  slip 

AM 

NE/SW 

9-26-96 

02:59:29 

2.8 

Strike  slip 

PM 

N/S 

11-1-96 

06:24:26 

2.7 

Strike  slip 

PM 

E/W 

12-20-96 

12:38:37 

2.6 

Dip  slip 

AM 

E/W 

2-14-97 

12:37:55 

2.5 

Dip  slip 

PM 

NW/SE 

3-14-97 

12:09:18 

3.6 

Dip  slip 

PM 

NW/SE 

5-25-97 

08:34:08 

2.6 

Dip  slip 

AM 

NW/SE 

06-27-97 

08:02:38 

3.5 

Strike  slip 

PM 

NW/SE 
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MODELING  OF  ROCK-BOLT  FLEXURAL  MODES 


Introduction 

Rock  bolts  are  important  components  in  modem  underground  mining  operations. 
Rock  bolts  provide  basic  support  for  unstable,  fractured  rock  by  coupling  rock  subunits  or 
linking  with  stable  rock  beyond  weak  rock  zones  (MESA,  1989).  Numerous  types  of 
rock  bolts  are  available  including  grouted  bolts,  mechanical  bolts,  and  friction  bolts.  The 
two  types  that  are  discussed  in  this  report  are  grouted  bolts  and  friction  bolts.  The 
Stillwater  Mine  in  Nye,  Montana  uses  grouted  bolts  that  are  eight  feet  long  and  split-set 
bolts  five  feet  long.  A  modal  analysis  of  the  split-set  bolt  is  discussed  in  a  following 
section.  Additionally,  computer  simulated  frequency  responses  of  a  ‘rebar’  bolt  are 
compared  with  closed-form  solutions  and  experimental  data. 

One  of  the  difficulties  in  determining  the  effectiveness  of  rock  bolts  is  the 
difficulty  of  determining  the  exact  condition  of  the  bolt  in  the  borehole  (Jumikis,  1983). 
One  goal  of  this  project  was  to  determine  how  close  computer-modeling  results  were  to 
closed-form  solutions  and  experimental  data  in  terms  of  frequency  response. 

Rock  Bolts 

An  overview  of  basic  rock-bolt  characteristics  is  contained  in  the  following 
sections. 

Split-set 

The  split-set  bolt  is  classified  as  a  friction  bolt  because  no  grouting  material  is 
used  to  hold  the  bolt  in  place.  This  friction  coupling  permits  the  bolt  to  yield  when 
movement  occurs  rather  than  fail  (Ingersoll-Rand,  1996).  The  bolt  is  a  hollow  tube, 
typically  1  'A”  in  diameter,  Vi”  thick,  with  a  A”  split  along  one  side  (Figure  25).  A  ring 
flange  is  welded  to  one  end  to  hold  the  faceplate  against  the  rock  and  the  other  end  is 
tapered  for  easier  installation. 

Boreholes  for  split-set  bolts  are  typically  1  3/8”  in  diameter.  The  bolt  is  driven 
into  the  borehole  with  a  driver  tool  and  instantly  exerts  a  radial  pressure  against  the  rock 
along  the  entire  length  of  the  bolt  (Ingersoll-Rand,  1996).  In  addition  to  the  radial 
pressure,  plate  loading  is  also  generated  immediately. 

The  natural  modes  of  the  split-set  bolt  were  calculated  using  finite  element 
analysis.  The  finite  element  mesh  that  was  used  was  a  rectangular  shell  mesh  with  a 
specified  thickness  of  0.10  inches.  The  mesh  was  five  feet  long,  which  is  the  length  of 
the  bolts  used  in  the  Stillwater  Mine.  A  close-up  of  the  mesh  is  shown  in  Figure  26. 
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Fig.  26.  Details  of  the  split-set  rock-bolt  mesh. 

In  Figure  26,  asterisks  represent  element  nodes.  The  rectangular  form  of  the 
elements  is  evident  and  the  boundary  conditions  used  in  the  modal  analysis  are  shown  as 
arrows  on  the  edge  of  the  bolt. 

Rebar 

The  rebar  bolt  represents  a  fundamentally  different  system  than  the  split-set  rock 
bolt.  For  the  case  of  rebar,  the  borehole  diameter  is  larger  than  the  rock-bolt  diameter. 
Typically  used  rock-bolt  diameters  are  Va”  and  common  borehole  diameters  are  1  3/8”. 
Grouting  material  choices  are  resin  or  cement.  Both  materials  have  decided  advantages 
and  disadvantages.  Resin  is  faster  setting,  but  is  more  expensive  than  cement. 

The  rock  bolt  is  installed  by  inserting  tubes  of  the  grouting  material  and  catalyst, 
then  driving  the  bolt  into  the  borehole.  Twisting  the  bolt  inside  the  borehole  breaks  the 
tubes  of  grouting  material  and  catalyst  and  mixes  them  together  to  form  an  adhesive  bond 
(Hoek,  1995).  Rebar  lacks  the  instantaneous  support  typified  by  split-set  bolts,  but  is 
nevertheless  an  effective  bolting  device. 

Because  a  rebar  bolt  is  solid;  unlike  a  split-set  bolt,  a  beam  mesh  was  used  in  the 
finite  element  method.  The  cross-sectional  area  of  the  bolt  varied  between  1 1/1 6”  and 
13/16”  because  of  the  ridges  present  on  rebar.  A  beam  mesh  was  created  by  creating  a 
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line  along  the  length  of  the  beam  and  specifying  a  cross-section.  The  element  length  was 
1/8”  which  is  equal  to  the  spacing  of  the  ridges. 

Split-set  Modal  Analysis 

One  of  the  primary  goals  of  the  modal  analysis  was  to  compare  computer¬ 
generated  results  with  in  situ  measurements  from  the  Stillwater  Mine.  Because  split-set 
bolts  were  not  used  in  the  mine  test  area,  the  split-set  analysis  was  performed  to  become 
familiar  with  the  modeling  software  and  to  compare  with  modeling  results  from  the  rebar. 
The  computer-generated  modes  are  found  in  Appendix  B  for  modes  one,  three,  five,  six, 
and  seven. 

The  split  along  the  side  of  the  hollow  tube  has  an  interesting  effect  on  the  modes 
of  the  split-set  bolt.  Mode  three  consists  of  a  radial  expansion,  along  vdth  a  slight 
bending  mode  (Figure  B-2,  Appendix  B).  The  first  mode  consists  of  a  simple  bending 
mode  with  little  or  no  radial  expansion  (Figure  B-1).  Modes  five  and  seven  display 
bending  modes  with  minimal  radial  expansion  (Figures  B-3,  B-5),  however  mode  six 
(Figure  B-4)  demonstrates  both  bending  and  expansion. 

Rebar  Bolt  Analysis 


Closed-form  solution  vs.  I-DEAS 

To  determine  the  validity  of  finite  element  modeling  of  rock  bolts,  the  results  of 
the  computer-simulated  solutions  were  compared  to  analytic  results.  Closed-form 
solutions  for  axial  vibrations,  lateral  vibrations,  and  a  clamped  string  were  used  for  the 
comparison. 


Axial  bolt  vibrations 

The  closed-form  solution  for  axial  vibrations  in  a  rock  bolt  that  is  clamped  or  free 
at  both  ends  is  given  by  (1)  (Tannant,  1995).  Calculations  from  this  equation  were 
compared  to  the  finite  element  modeling  results  for  the  same  bolt  configuration.  In 
equation  (1),  f„  is  the  axial  frequency  of  the  rock  bolt,  Y  is  Young’s  modulus  of  steel,  p  is 
the  density  of  steel,  and  L  is  the  length  of  the  rock  bolt. 


V  p 

21  ; 


(1) 


Young’s  modulus  for  steel  is  200  GPa  and  the  density  of  steel  is  7800  kg/m^  Using  a 
length  of  two  meters  for  the  bolt  yields  a  first  harmonic  frequency  of  1270  Hz. 

Using  I-DEAS,  a  beam  mesh  two  meters  long  was  created  with  a  diameter  of  16 
mm.  Using  fixed  boundary  conditions  at  both  ends,  the  first  axial  mode  of  the  bolt  was 
calculated.  The  frequency  of  the  first  harmonic  frequency  of  the  bolt,  as  calculated  by  I- 
DEAS,  was  1290  Hz.  The  percent  difference  between  the  I-DEAS  solution  and  the 
closed-form  solution  is  1 .57%.  The  two  solutions  are  very  close  indicating  that  I-DEAS 
accurately  models  the  axial  frequency  of  a  bolt  clamped  at  both  ends. 
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Lateral  bolt  vibrations 

The  second  closed-form  equation  is  for  lateral  bolt  vibrations  and  is  given  in  (2) 
where  f„  is  the  fundamental  frequency,  L  is  the  length  of  the  rock  bolt,  Y  is  Young’s 
modulus  for  steel  (200  GPa),  p  is  the  density  of  steel  (7800  kg/m3),  I  is  moment  of 
inertia,  and  A  is  the  circular  cross-sectional  area  of  the  bolt. 


As  for  (1),  the  length  of  the  rock  bolt  is  two  meters  and  the  diameter  is  16  mm.  Moment 
of  inertia  (I)  is  calculated  using  (3)  (Weaver,  Jr.,  1980). 


1  = 


4 


(3) 


The  variable  is  a  boundary  condition  constant  the  values  for  which  were  found  in 
Engineering  Vibration  (Inman,  1996).  Values  for  each  boundary  condition  and  the 
fundamental  frequencies  from  both  solutions  are  given  in  Table  3. 


Table  3.  Lateral  Bolt  Vibration  Fundamental  Frequencies 


Boundary 

Condition 

Closed-form  solution 
(Hz) 

I-DEAS  solution 
(Hz) 

% 

difference 

Fixed-Free 

3.52 

2.84 

2.88 

1.41 

Fixed-Hinged 

15.4 

12.41 

12.62 

1.69 

Fixed-Fixed 

22 

17.73 

18.31 

3.27 

The  percent  differences  show  that  the  I-DEAS  finite  element  modeling  accurately  models 
the  lateral  bolt  vibrations. 

Clamped  string 

The  third  closed-form  equation  was  for  the  fundamental  frequency  of  a  string 
clamped  at  both  ends  under  tension  (4)  (Tannant,  1995)  where  is  the  frequency,  L  is  the 
length  of  the  rock  bolt  (two  m),  T  is  the  tension  in  the  string,  p  is  the  density  of  steel 
(7800  kg/m3),  and  A  is  the  cross-sectional  area  of  the  rock  bolt. 

The  diameter  of  the  ‘string’  is  16  mm.  Typical  tensions  in  rock  bolts  range  from  10  kN  to 
50  kN.  In  I-DEAS,  the  minimum  and  maximum  of  the  range  were  used  to  calculate 
frequencies.  The  results  are  given  in  Table  4. 
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Table  4.  Tensions  and  Frequencies  for  Closed-Form  Solution  Comparison  with  Finite 

Element  Modeling 


Tension 

Closed-form  solution 
(Hz) 

I-DEAS  solution 
(Hz) 

%  difference 

10  kN 

20 

29 

45.00 

50  kN 

45 

51 

1  n  o  n 

1  J.J 

The  figures  showing  the  results  above  are  given  in  Appendix  C  as  C-1  and  C-2.  The 
percent  differences  from  this  comparison  are  high  and  can  be  explained  by  the  following: 
the  equation  for  the  clamped  string  was  for  a  non-rigid  string.  However,  a  rock  bolt  is 
rigid  and  this  most  likely  accounts  for  the  large  difference  values. 

Results 

The  results  from  the  axial  bolt  vibration  and  lateral  bolt  vibration  were 
encouraging.  Percent  differences  for  the  first  two  comparisons  were  low  and  indicative 
of  realistic  modeling  parameters.  The  percent  differences  for  the  clamped  string  were 
large  but  can  be  explained  by  the  difference  in  rigidity  of  the  rock  bolt  and  the  string 
closed  form  solution. 

Experimental  data  vs.  I-DEAS  modeling 

Another  experiment  was  conducted  to  determine  how  accurately  I-DEAS  models 
compare  with  experimental  results. 

Laboratory  testing 

The  experiment  consisted  of  putting  two  rock  bolts  under  tension  and  tapping  the 
ends  of  the  rock  bolts  with  a  hammer.  Using  a  microphone,  the  frequencies  of  the 
induced  vibrations  were  recorded.  The  experiment  was  semi-controlled  in  the  sense  that 
the  tensions  were  measured,  but  not  the  amount  of  force  for  each  tap. 

The  applied  tensions  to  both  the  rebar  and  cold-rolled  steel  bolts  were  1000 
pounds  and  5000  pounds.  Each  bolt  was  tapped  twice  on  the  top  and  the  bottom.  The 
sound  recordings  were  transformed  to  the  frequency  domain  and  plotted.  Signal-to-noise 
ratios  were  also  calculated. 

I-DEAS  modeling 
Cold-rolled  steel 

The  initial  and  most  important  part  of  modeling  in  I-DEAS  is  to  create  the  correct 
finite  element  mesh.  Grip  marks  were  imbedded  into  the  bolt  from  the  experiment,  and 
the  length  of  each  grip  was  3  5/8”  long.  The  top  grip  started  one  inch  from  the  top  of  the 
bolt,  and  the  bottom  grip  was  7  'A”  from  the  base  of  the  bolt.  Total  length  of  the  rock 
bolt,  from  the  extra  mass  around  the  eye,  was  19  VC,  leaving  an  effective  length  between 
the  grips  of  4”.  The  diameter  of  the  rock  bolt  was  5/8”  (Figure  27). 
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A  beam  cross-section  of  5/8”  diameter  was  created  for  the  entire  beam.  The 
elements  were  formed  to  be  1/8”  long.  To  recreate  the  effects  of  the  grips,  boundary 
conditions  were  set  at  every  node  in  the  areas  that  the  grips  were  placed.  All  of  the 
rotational  degrees  of  freedom  and  translation  in  the  axial  direction  were  free  in  the 
boundary  conditions  (Figure  27). 


Fig.  27.  Cold-rolled  steel  bolt  mesh. 

In  addition  to  the  boundary  conditions,  a  load  of  1000  pounds  was  applied  to  the 
boundary  conditions  closest  to  each  end  to  represent  the  tension  applied  to  the  rock  bolt. 
To  compare  results  from  the  sound  recordings  and  modeling  in  I-DEAS,  solutions  were 
developed  for  the  normal  modes  of  the  bolt.  These  normal  modes  represent  what  occurs 
if  the  bolt  is  excited  at  a  certain  frequency  (Figure  28). 


Fig.  28.  Normal  mode  at  4240  Hz  for  cold-rolled  steel  bolt. 

The  normal  mode  frequency  for  the  cold-rolled  steel  bolt  under  a  tension  of  1000 
pounds  was  4240  Hz  (Figure  28).  This  agrees  with  the  experimental  results  shown  later. 
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Rebar 


A  similar  analysis  was  completed  for  a  rebar  bolt.  The  dimensions  of  the  rebar 
bolt  were  different  than  the  cold-rolled  steel  bolt.  The  total  length  of  the  rebar  was  26 
5/8”.  The  grips  were  3  5/8”  from  one  end,  and  2  7/8”  from  the  other  end.  As  with  the 
cold-rolled  steel,  the  grip  lengths  are  3  5/8”  but  they  react  in  different  ways.  Because  the 
rebar  has  ridges,  the  grip  contact  points  occur  only  along  the  ridges.  Originally,  the  bolt 
model  was  meshed  by  making  the  element  widths  the  same  width  as  the  ridge  spacing 
(1/8”)  with  a  cross-section  1 1/16”  in  diameter.  However,  this  was  changed  so  that  the 
diameter  of  every  third  element  was  13/16”,  which  is  the  diameter  of  the  ridges  (Figure 
29). 


Fig.  29.  Rebar  rock  bold  mesh  for  finite  element  modeling. 

The  boundary  conditions  were  applied  to  the  bolt  along  the  nodes  of  the  ridge 
elements.  Again,  all  but  the  end  boundary  condition  allowed  free  rotational  displacement 
and  free  translational  displacement  in  the  axial  direction.  The  end  boundary  condition 
allowed  free  rotational  displacement,  but  all  of  the  translational  displacements  were  fixed 
(Figure  30). 


Fig.  30.  Solid  rebar  and  boundary  conditions. 

Next,  an  excitation  function  was  created  in  I-DEAS  to  provide  an  impulse  to  the 
bolt  to  simulate  a  hammer  tap  on  the  experimental  bolt.  The  applied  excitation  was  a  unit 
force  of  one  pound.  Because  the  actual  force  of  the  hanuner  tap  could  not  be  easily 
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determined,  a  unit  force  of  one  pound  was  used  for  the  excitation.  The  force  was  applied 
axially  and  the  reactions  that  are  plotted  are  frequency  responses  in  the  axial  direction. 
Without  completely  controlled  experimental  conditions,  the  hits  may  have  been  partially 
lateral,  resulting  in  bending  modes  that  do  not  appear  on  the  computer  simulated  response 
plot.  Results  of  the  excitation  are  shown  in  Figure  3 1 . 
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Fig.  31.  Frequency  response  —  rebar  rock  bolt  (bottom  excitation). 

The  vertical  axis  of  the  graph  represents  the  magnitude  of  the  displacement  even 
though  I-DEAS  uses  the  term  modulus  as  a  general  term  for  all  vertical  axis  values. 
Below  the  horizontal  -axis  is  a  list  of  the  rock-bolt  nodes  at  which  each  displacement  was 
calculated.  When  the  bolt  was  struck  in  the  lab,  the  sound  that  was  recorded  is  a  result  of 
all  of  the  nodes.  In  I-DEAS,  the  contribution  from  all  the  nodes  is  estimated  by  choosing 
a  small  number  of  nodes  along  the  length  of  the  bolt.  The  first  axial  mode  in  Figure  3 1  is 
shown  as  a  peak  at  21 50  Hz.  Experimental  results  correlate  with  the  computer  analysis 
for  a  rebar  bolt  excited  on  the  bottom.  We  used  the  same  excitation  function  on  the  bolt 
top.  The  results  were  shifted  slightly  from  2150  Hz  to  2200  Hz  from  the  bottom 
excitation  results  (Figure  32).  This  is  expected  because  the  length  of  the  bolt  does  not 
change,  nor  do  the  boundary  conditions  on  the  bolt.  The  only  thing  that  changes  is  the 
length  from  each  end  to  the  grips,  but  that  change  is  small  (3/4”). 
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Fig.  32.  Frequency  response  -  rebar  rock  bolt  (top  excitation). 

Results 

The  results  from  the  comparison  of  experiment  results  and  computer  modeling  are 
encouraging.  Computer  modeling  is  shown  to  be  reliable  when  analyzing  the  effects  of 
excitations  and  forces  on  rock  bolts.  The  same  initial  axial  modes  that  are  appearing  in 
the  experimental  recordings  appear  in  the  computer  produced  frequency  response. 
Because  the  computer  modeling  compares  well  to  experimental  and  closed-form 
solutions,  it  will  next  be  applied  to  the  rebar  bolts  found  in  the  Stillwater  Mine. 

Stillwater  Mine  Experiments 

To  simulate  the  rock-bolt  conditions  in  the  Stillwater  Mine,  a  finite  element  mesh 
was  constructed  that  was  eight  foot  long.  Two  varying  sets  of  boundary  conditions  were 
used  with  constant  excitations  for  comparison  with  experimental  results. 

The  eight-foot  rebar  was  meshed  the  same  way  that  the  smaller  rebar  bolt  was 
meshed  in  the  prior  simulation  using  the  ridges  along  the  bolt.  In  a  previous  report  by  the 
U.S.  Bureau  of  Mines,  a  finite  element  mesh  was  used  to  model  a  rock  bolt,  however, 
those  results  did  not  agree  with  experimental  results  (Serbousek,  1987).  One  of  the 
reasons  given  for  the  disparity  was  that  the  finite  element  model  did  not  truly  match  the 
behavior  of  the  rock  bolt  because  the  ridges  of  the  bolt  were  not  modeled.  In  order  to 
avoid  this  problem,  ridges  were  included  in  our  modeling.  Boundary  conditions  were 
used  at  each  node,  except  for  the  last  four  inches,  to  represent  a  completely  bonded  bolt. 
In  the  Stillwater  Mine,  part  of  the  rebar  protrudes  from  the  wall,  which  is  why  four  inches 
of  the  bolt  are  not  restrained. 

Meshing  the  bar  and  applying  the  boundary  conditions  provides  a  bolt  that  is 
completely  bonded,  even  in  between  the  ridges.  In  the  modeled  rebar  mesh,  .the  boundary 
conditions  were  along  the  ridges.  This  is  different  from  the  mine  conditions  in  that  the 
Stillwater  Mine  bolts  are  grouted,  so  that  the  grout  bonds  to  all  of  the  rock-bolt  surface, 
not  only  to  the  ridges.  The  results  from  an  impulse  excitation  are  shown  in  Figure  33. 
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Fig.  33.  Modeled  Stillwater  Mine  rebar  rock-bolt  frequency  response  - 
completely  bonded. 

Frequency  transformed  sound  recordings  from  the  Stillwater  Mine  rebar  rock  bolts 
exhibit  a  similar  series  of  peaks;  however,  the  amplitude  of  the  peaks  decreases  with 
frequency  and  the  location  of  the  peaks  is  shifted.  The  variations  in  the  mine-recorded 
frequencies  from  the  modeled  frequencies  could  result  from  gaps  in  bonding,  excitation 
of  lateral  as  well  as  axial  modes,  or  non-impulse  type  excitation. 

A  second  set  of  boundary  conditions  was  used  for  a  rebar  rock  bolt  that  was  half- 
bonded.  The  completely  bonded  end  of  the  bolt  represents  the  portion  of  a  rock-bold 
inside  a  rock  mass  and  the  unbonded  end  represents  the  portion  of  the  rock-bold  external 
to  the  rock.  The  results  from  the  frequency  modeling  analysis  for  the  half-bonded  bolt 
are  shown  in  Figure  34. 


Fig.  34.  Modeled  Stillwater  Mine  rebar  rock-bolt  frequency  response  -  half- 
bonded. 
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The  most  noticeable  difference  between  the  completely  bonded  bolt  results  and 
the  half-bonded  bolt  results  is  the  difference  in  the  number  of  peaks.  This  seems  to 
indicate  that  the  degree  of  bonding  is  related  to  the  number  of  peaks  in  the  frequency 
response  of  the  bolts.  As  can  be  seen,  the  first  peak  of  the  completely  bonded  bolt  is  at 
the  same  frequency  as  the  peak  of  the  half-bonded  bolt.  Further  experiments  and  research 
can  identify  bonding  conditions  of  rock  bolts  in  more  detail.  One  possibility  is  using  a 
pulse-echo  system  to  measure  the  travel  times  of  sonic  pulses  in  bolts  (Tadolini,  1990) 
rather  than  a  microphone. 


Conclusions 

The  purpose  of  this  project  was  to  determine  the  feasibility  of  modeling  rock  bolts 
and  their  associated  frequency  responses  with  finite  element  modeling  as  found  in  the 
computer  software  I-DEAS.  When  compared  to  closed-form  solutions,  1-DEAS  modeled 
axial  bolt  vibrations  and  lateral  bolt  vibrations  accurately  (within  4%  for  all  cases). 

Laboratory  experimental  results  were  also  compared  to  computer  modeled 
frequency  responses.  Again,  the  results  of  laboratory  experiments  and  the  computer 
models  were  similar.  This  modeling  was  also  applied  to  the  rebar  bolts  used  in  the 
Stillwater  Mine.  Two  situations  were  modeled:  1)  completely  bonded  bolts,  and  2)  half- 
bonded  bolts.  The  results  of  both  computer-simulated  situations  were  compared  to  actual 
recordings  obtained  from  the  Stillwater  Mine.  It  can  be  inferred  that  the  sample  bolts  in 
the  Stillwater  Mine  were  well  bonded. 

In  conclusion,  the  modeling  of  rock  bolts  using  the  finite  element  methods  as  in 
the  computer  software  I-DEAS  is  a  useful  technique.  Further  research  will  lead  to  a  more 
detailed  understanding  of  rock-bolt  bonding  and  its  role  in  mining  safety. 
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FREQUENCY  RESPONSE  OF  ROCK-BOLT  SOUNDINGS  FROM 
HAMMER  STRIKES 

Introduction 

This  research  consisted  of  analysis  of  acoustic  rock-bolt  signatures  to  determine 
frequency  content.  Rock  bolts  relieve  loading  stress  and  add  support  along  mine  walls.  Under 
high  stress  loading,  rock  bolts  can  loosen  and  fail.  This  can  trigger  cave-ins  during  blasting, 
earthquakes,  and  rockbursts  (Otuonye,  1998).  If  viable,  this  research  would  contribute  to  an 
inexpensive,  quick  solution  to  locate  areas  of  increased  loading  and  possible  hazard. 

The  idea  for  the  experiment  stems  from  miner’s  folklore.  Underground  miners 
maintain  that  they  can  predict  mine  wall  loading  by  the  sound  of  a  rock  hammer  hitting  a  rock 
bolt.  To  investigate  the  merit  of  this  belief,  recordings  of  rock  hammers  striking  rock  bolts 
were  made  on  a  pillar  under  investigation  in  the  Stillwater  Mine.  To  provide  a  better 
understanding  of  the  variables  involved,  rock  bolts  were  also  analyzed  in  a  laboratory 
environment  under  controlled  amounts  of  loading. 

Experimental  Procedure 

A  two-channel  Sony  ECM-929L  microphone  was  used  for  the  acoustic  data  collection 
for  this  project.  The  properties  of  this  microphone  include  a  dynamic  range  of  more  than  94 
dB  and  signal-to-noise  ratio  greater  than  48  dB. 

Laboratory  experiment 

To  detect  changes  in  frequency  content,  experiments  were  performed  on  the  Montana 
Tech  campus  March  13, 1998.  These  experiments  were  performed  on  a  DYWIDAG  rebar 
rock  bolt,  donated  by  the  Stillwater  Mine,  and  a  Simmons  eye  bolt,  donated  by  J.P.  Simmons 
Iron  Works  Inc.  (Figure  35). 

The  DYWIDAG  rebar  rock  bolt  is  composed  of  %  inch  #6  grade  70  steel  with  1 3/1 6- 
inch  ridges.  This  rock  bolt  is  identical  to  those  used  in  the  Stillwater  Mine.  The  Simmons  eye 
bolt  is  constructed  of  5/8  inch  cold-rolled  steel.  The  eight-foot  DYWIDAG  bolt  was  cut  to  a 
length  of  27  inches,  approximately  the  same  length  as  the  Simmons  eye  bolt.  Each  of  these 
Rock-bolts  was  then  inserted  into  the  Tinius-Ohlsen  tension/compression  tester.  The 
Simmons  eye  bolt  was  secured  by  clamps  directly  above  the  eye  and  directly  above  the  bolt 
threads.  The  DYWIDAG  rock  bolt  was  secured  approximately  3  inches  from  either  end 
(Figure  36). 
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Fig.  36.  Rock-bolt  and  machine  setup. 


The  tensional  breaking  point  of  each  rock  bolt  was  then  calculated.  This  helped  to 
determine  the  range  of  tensions  under  which  each  rock  bolt  would  be  placed.  The  results  of 
calculations  are  as  follows: 

Simmons  eye  bolt  (Oberg  et  al.,  1996): 

tensional  strength  =  60,000  xn 

in 

=  18,408  Ih 

DYWIDAG  rock  bolt  (Oberg  et  al.,  1996): 

Ih 

tensional  strength  =  20,000  xtt 

in' 

=  8,836  lb 

Using  the  tensional  strength  as  the  maximum  tension  the  rock  bolt  can  withstand, 
measurements  were  taken  for  the  DYWIDAG  rock  bolt  at  1000, 2000, 3000, 4000,  and  5000 
lb.  Measurements  for  the  Simmons  eye  bolt  were  taken  at  tensions  of  1000, 2000, 3000, 

4000,  5000,  10000, 15000  lb.  The  tension  ranges  lie  within  the  approximately  2000  to  1 1000 
lb.  tensions  found  in  a  typical  mine  site  (Means,  1976). 

Stillwater  Mine  data  acquisition 

The  data  for  this  project  were  acquired  at  the  Stillwater  Mine.  Acoustic  data  were 
recorded  at  the  Bin  #2  Pillar  at  the  3200-foot  level.  The  station  locations  coincide  with  those 
used  for  the  ground  penetrating  radar  tomographic  study.  The  station  locations  and  pillar 
geometry  are  shown  in  Figure  1 .  Each  station  contains  an  upper  and  lower  rock  bolt.  Each  of 
these  individual  rock  bolts  was  struck  on  the  end  and  side  two  times.  The  sound  produced  by 
each  rock  bolt  was  recorded  and  subsequently  analyzed  in  the  frequency  domain. 

Expected  Outcome 

Some  uncertainty  concerning  the  true  Young’s  modulus  of  rebar  arose  from  a  paper 
by  Tadolini  (1990).  His  research  presented  some  useful  information  for  calculating  Young’s 
modulus  for  rebar.  Tadolini  found  that  a  5/8-inch  diameter,  4-foot  long  rebar  rock  bolt 
stretches  0.0064  inches  per  100  lb.  Using  this  information.  Young’s  modulus  was 
recalculated  to  present  a  range  that  might  more  accurately  represent  rebar.  For  ease  of 
calculations,  English  units  were  converted  here  to  metric  units. 

tensile  strain  =  =  1.333  x  10"^ 

tensile  stress  =  ^  =  2.24  x  10* 

A 

Young's  modulus  =  stress  =  1.68  x  lO"  Pa 

°  tensile  strain 

=  2.43  X  10*-% 
in 
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To  verify  the  Young’s  modulus  recalculations  for  rebar,  an  experiment  was  conducted 
to  experimentally  determine  Young’s  modulxis  for  #6  grade  70  steel  (rebar).  A  rod  of  V2  inch 
diameter  rebar  was  placed  under  increasing  tension  in  the  Tinius-Ohlsen  testing  machine.  The 
change  in  length  was  measured  with  an  extensiometer.  The  results  are  presented  in  Table  5. 

Table  5.  Young’s  Modulus  Laboratoiy  Measurements  for  Rebar 


aL  (X 

Slope 

Load  (lb) 

0.001  in)  Force/Area 

aL/L 

(lb/in'‘2) 

0 

0.00 

0 

0 

O.OOOOE+00 

1000 

1.50 

5093 

0.0015 

3.3953E+06 

2000 

3.50 

10186 

0.0035 

2.9103E+06 

3000 

5.50 

15279 

0.0055 

2.7780E+06 

4000 

7.20 

20372 

0.0072 

2.8294E+06 

5000 

9.50 

25465 

0.0095 

2.6805E+06 

6000 

11.10 

30558 

0.0111 

2.7530E+06 

7000 

13.00 

35651 

0.013 

2.7424E+06 

8000 

15.00 

40744 

0.015 

2.7162E+06 

9000 

17.00 

45837 

0.017 

2.6963E+06 

10000 

19.00 

50930 

0.019 

2.6805E+06 

Breaking 

Young's 

strength 

modulus 

= 

15880  lbs 

= 

2.5620E+06 

Using  linear  regression,  a  best-fit  slope  was  calculated  from  the  experimental  data  to 
determine  Young’s  modulus.  This  value  was  approximately  2.56  x  lOMb  in'".  For 
comparison,  Young’s  modulus  for  cold  rolled  steel  is  2.9  x  10^  lb  in'^  (Sears  et  al.,  1987).  This 
calculated  value  of  Young’s  modulus  for  #6  Grade  70  steel  was  used  in  all  following 
calculations  pertaining  to  the  DYWIDAG  (rebar)  frequencies. 

By  striking  a  rock  bolt  on  both  the  side  and  end,  two  types  of  fundamental  frequencies 
are  produced.  Striking  the  end  activates  axial  (longitudinal)  modes,  while  striking  the  side 
activates  bending  (transverse)  frequencies  (Morse,  1948).  Tannant  et  al.  (1995)  report  on 
experiments  similar  to  this  investigation  on  rock  bolts.  In  his  experiments,  strain  gauges  were 
placed  on  either  side  of  the  bolt.  When  struck,  one  side  of  the  bolt  was  in  compression  while 
the  other  was  in  tension.  Based  on  this  concept,  Tannant  et  al.  viewed  the  rock  bolt  as  an 
anchored  beam.  Had  the  rock  bold  exhibited  uniform  deformation,  the  rock  bolt  would  more 
closely  approximate  an  anchored  string  (Morse  1948).  The  calculations  for  the  approximate 
fundamental  frequencies  are  as  follows: 
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Axial  fundamental  frequencies  (Morse,  1948): 


where 


Fp  =  compressional  velocity  ■ 

Y  =  Young’s  modulus  for  steel, 

a 

p  =  density  of  steel  =  7.8  — ?, 

^  cm^ 

n  =  frequency  mode  number,  and 
L  =  length  of  the  rebar. 


Y 
P  ’ 


DYWIDAG  fundamental  frequency  =  2745  Hz 

Simmons  eye-bolt  fundamental  frequency  =  3940  Hz 

Stillwater  Mine  rock-bolt  (rebar)  fundamental  frequency  =  1 26 1  Hz 


Bending  fundamental  frequencies  (Morse,  1948): 


IL 

Ap 


where 

nd'’ 

I  =  moment  of  inertia  =  , 

o4 

Y  =  Young’s  modulus  for  steel, 

a 

p  =  density  of  steel  =  7.8 

A  =  cross  sectional  area, 

L  =  length  of  the  rebar, 

A„  =  3 .52  for  one  fixed  end  and  one  free  end, 

1 5.4  for  one  fixed  end  and  one  hinged  end, 

22.0  for  both  ends  fixed,  and 
d  =  diameter  of  the  rebar. 

DYWIDAG  fundamental  frequency  =  10  to  84  Hz 

Simmons  eye-bolt  fundamental  frequency  =  15  to  94  Hz 

Stillwater  Mine  rock-bolt  (rebar)  fundamental  frequency  =  4  to  28  Hz 

Mine  Recording  Data  Analysis 

All  of  the  microphone  recordings  were  converted  to  digital  form.  A  high  sampling 
rate  of 44100  Hz  sampling  rate  was  used  to  prevent  aliasing  and  preserve  high  frequencies. 
Using  the  Sound  Editor  application  on  a  Silicon  Graphics  Inc.  (SGI)  workstation,  each 
individual  rock-bolt  hit  was  extracted  as  a  digital  data  subset  from  the  analog  recording. 
These  sound  files  were  then  converted  from  aij^  files  to  au  files. 
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A  Matlab  program  was  developed  to  analyze  the  sound  files  in  the  frequency  domain. 
The  program  begins  by  separating  the  recordings  of  the  two  channels.  Next,  the  user 
graphically  picks  noise  and  signal  zones  (Figure  37).  / 


Fig.  37.  Example  of  single  hammer  hit  recording. 

Finally,  the  power  spectral  content  of  both  the  signal  and  noise  intervals  were  calculated. 
These  values  were  used  to  calculate  signal-to-noise  ratio  (S/N).  Each  hammer-hit  recording 
was  analyzed  in  a  similar  fashion. 

S/N  averaged  approximately  46  dB  for  all  of  the  recordings  in  the  higher  frequency 
ranges.  S/N  was  significantly  lower  in  the  lower  frequency  range.  This  is  unfortunately  the 
zone  where  the  bending  fundamental  frequencies  are  located  (Figure  38).  For  this  reason, 
research  focused  upon  axial  measurements. 

Lab  experimental  data 

Two  frequency  peaks  at  approximately  500  Hz  and  750  Hz  were  prevalent  on  all  of 
the  laboratory  recordings.  To  expleiin  these  frequencies,  a  plot  of  the  Tinius-Ohlsen  machine 
noise  was  made  (Figure  39).  This  plot  clearly  demonstrates  that  the  two  frequency  peaks  are 
the  result  of  machine  noise.  The  Simmons  eye-bolt,  axial,  bottom-hit  results  are  shown  in 
Figure  40. 
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Power  (dB)  Signal  to  Noise  Ratio 


Fig.  38.  Typical  low  signal-to-noise  ratio  (S/N)  of  lower  frequencies. 


Fig.  39.  Tinius-Ohlsen  machine  noise  frequencies. 


53 


Fig.  40.  Axial,  bottom-hit  Simmons  eye-bolt  results. 

The  calculated  fundamental  frequency  of  the  rock  bolt  including  the  eye  is  3503  Hz. 
This  value  did  not  correspond  with  the  approximately  4200  Hz  peak  observed  in  the 
laboratory  data.  However,  using  this  higher  frequency  value,  a  calculated  rock-bolt  length  of 
22.5  inches  is  obtained.  This  length  corresponds  to  the  rock  bolt  without  the  eye. 

The  peak  located  at  approximately  1 500  Hz  is  theorized  to  be  a  combination  of  axial 
and  bending  frequencies  into  a  torsional  frequency  (Maxwell,  1941).  Both  a  non-direct  end 
hit  and  placement  of  the  Tinius-Ohlsen  clips  could  cause  this  result.  The  clips  contain  notched 
teeth,  and  were  secured  approximately  three  inches  from  each  bolt  end.  Previously  work 
shown  earlier  in  this  report  on  rock-bolt  finite  element  modeling  shows  computer  models  that 
support  this  premise.  Both  axial  and  bending  components  were  included  in  the  computer 
simulation,  and  frequency  peaks  similar  to  the  1 500  Hz  peak  appear  in  the  results. 

The  axial  top  Simmons  eye-bolt  recordings  display  two  prominent  peaks  (Figure  41). 
The  peak  at  3450  Hz  was  close  to  the  calculated  value  and  when  the  length  was  calculated 
using  the  observed  frequency  value,  a  length  of  28  inches  was  obtained.  This  length 
corresponds  to  the  length  of  the  rock  bolt  including  the  eye.  The  peak  at  approximately  2250 
Hz  is  again  theorized  to  be  a  combination  of  the  axial  and  bending  frequencies.  The  piece  of 
DYWIDAG  rebar  rock  bolt  from  the  Stillwater  Mine  was  placed  under  varying  amounts  of 
tension  and  the  frequency  content  was  determined.  However,  none  of  the  results  from  this 
analysis  exhibited  the  calculated  fundamental  frequency. 

The  frequency  analysis  results  of  the  axial  bottom  DYWIDAG  results  are  shown  in 
Figure  42  and  the  results  from  the  axial  top  DYWIDAG  rebar  are  shown  in  Figure  43. 
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Frequencies  0  to  5000  Hz 
Top  Axial  DYWIDAG  Results 


Frequency  (Hz) 

Fig.  43.  Axial  top  DYWIDAG  rebar  rock-bolt  results. 

Two  dominant  frequency  peaks  occur  at  approximately  1250  Hz,  and  750  Hz.  Neither  of 
these  frequencies  is  near  the  calculated  fundamental  frequency.  To  explain  these  results,  a 
finite  element  computer  model  was  based  upon  this  data  set.  Incorporated  into  the  model  are 
the  13/16-inch  ridges  along  the  bolt,  the  position  of  the  clamps,  and  boundary  conditions  of 
the  bolt  upon  these  ridges.  An  axial  impulse  was  simulated  on  this  model.  Modeled  results 
give  a  peak  at  approximately  2100  Hz.  Frequencies  of  approximately  1250  Hz  and  750  Hz 
are  also  excited  when  stuck  by  an  axial  hit.  The  combination  of  clamping  upon  the  ridges  and 
the  ridges  themselves  has  reproduced  these  lower  frequencies  in  the  computer  model. 

Stillwater  Mine  data  analysis 

The  fundamental  frequency  expected  in  the  Stillwater  Mine  rock-bolt  data  set  was 
calculated  to  be  1260  Hz.  This  frequency  appeared  on  almost  every  axial  sound  recording 
from  the  mine  (Figure  44).  However,  at  some  of  the  stations,  the  dominant  frequency  was 
located  at  approximately  400  Hz.  This  large  variation  from  the  calculated  fundamental 
frequency  is  difficult  to  explain.  The  difference  might  be  due  to  different  length  rock  bolts, 
variable  coupling  of  the  bolt  to  the  rock,  or  non-direct  hits  on  the  ends  of  the  bolts.  In 
addition,  the  plates  securing  the  rock  bolts  to  the  rock  face  exhibited  a  large  degree  of 
variation  in  tightness.  Also  note  that  the  finite  element  model  of  this  rock  bolt  produced 
frequencies  at  approximately  2200  Hz  and  3300  Hz.  These  frequencies  are  also  observed  in 
the  mine  recordings  (Figure  44). 


56 


Stillwater  Mine  Axial  Data 


Fig.  44.  Stillwater  Mine  axial  sound  recordings  analysis  results. 

A  set  of  tables  allows  easy  comparison  of  frequency  peaks  contained  in  the  various 
data  sets.  Table  6  lists  the  frequency  peaks  for  the  laboratory  experiments  for  the  DYWIDAG 
and  Simmons  eye  rock  bolts.  Table  7  lists  several  peak  frequencies  for  typical  sound 
recordings  in  the  Stillwater  Mine. 

Table  6.  Dominant  Frequencies  for  Laboratory  Measurements 


1000  Ib  tension 

5000/15000  lb  tension 

bottom  axial  Simmons  eye 

4157  Hz  and  1530  Hz 

4226  Hz  and  1578  Hz 

top  axial  Simmons  eye 

3487  Hz  and  2300  Hz 

3491  Hz  and  2412  Hz 

bottom  axial  DYWIDAG 

2070  Hz 

2078  Hz 

top  axial  DYWIDAG 

1187  Hz  and  750  Hz 

1 188  Hz  and  765  Hz 
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Table  7.  Dominant  Frequencies  for  Stillwater  Mine  Measurements 


Station 

Frequency  (Hz) 

Two 

1362 

Four 

1230 

Five 

1227 

Seven 

1250 

Ten 

1412 

Fourteen 

1311 

Conclusions 

The  bending  or  longitudinal  frequencies  for  this  data  set  are  very  low.  As  such,  they 
peak  in  an  area  with  a  low  signal  to  noise  ratio.  Therefore,  the  focus  of  the  research  was  on 
the  axial  or  transverse  frequencies.  The  background  noise  produced  by  the  Tinius-Ohlsen 
machine  was  plotted.  Frequencies  of  approximately  550  Hz  and  750  Hz  were  produced  by 
this  machine. 

Using  the  Simmons  eye  rock  bolt,  frequencies  representing  the  calculated  fundamental 
frequency  were  apparent  in  the  data.  When  striking  the  bottom  of  this  bolt,  it  is  postulated  that 
the  wave  travels  up  the  bolt  to  the  point  where  the  attachment  ring  begins  and  then  returns. 

The  approximate  4200  Hz  frequency  peak  corresponds  with  this  result.  Hitting  the  top  of  the 
rock  bolt  upon  the  attachment  ring,  the  resulting  wave  would  have  to  travel  a  longer  distance. 
This  distance  of  approximately  28  inches  corresponds  with  the  3200  Hz  observed  frequency. 
Also,  at  2250  Hz,  frequency  peaks  were  visualized.  Our  theory  about  these  peaks  concerns 
the  result  of  both  axial  and  bending  frequencies  effecting  the  bolt.  These  could  be  due  to  the 
placement  of  the  clamps  on  the  Tinius-Ohlsen  machine  and/or  non  direct  hits  on  the  end  of  the 
bolt.  Our  theory  is  supported  by  the  computer  models. 

Two  dominant  frequencies  were  evident  in  the  DYWIDAG  rebar  results  at  750  Hz 
and  1250  Hz.  These  frequency  peaks  do  not  correspond  to  any  calculated  fundamental 
frequencies.  However,  computer-modeling  results  suggest  possible  explanations.  The 
frequency  peaks  could  result  from  the  ridges  on  the  rock  bolt,  the  securing  of  the  rock  bolt  off 
the  ends  of  the  bolt,  and  the  securing  of  the  rock  bolt  upon  the  ridges. 

The  Stillwater  Mine  data  exhibited  the  fundamental  frequency  peak  at  approximately 
1260  Hz  over  a  range  of 400  Hz.  This  somewhat  wide  range  could  result  from  non-direct 
strikes  on  the  ends  of  the  rock  bolts,  different  rock-bolt  lengths,  poor  coupling  of  the  bolt  to 
the  rock,  or  the  ridges  on  the  sides  of  the  DYWIDAG  rebar  rock  bolts. 
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Appendix  A 


Rockburst  Events 
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185  45  -90  .00  12  8.05  .72  .00  13  8  25 
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95091905 
215  70  lliO 


UFDWN 

00  DILflTflTION 
+  +  COMPRE9SI0N 


.0047-26.10115-43.09-20.10  .30  .0 

.08  12  8.05  .60  .00  13  5  10 
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95091903 
150  50  -80 


.30 


.0 


.00  47-26.10  115-43.09  -20.10 
.00  10  8.05  .63  .00  3  0  0 
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.0 


96020800  .00  117-26.10  IIS-US. 09  -22.22  .60 

2SQ  SQ-IHQ  .09  11  9.05  .69  .00  20  5  lOCiK 
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96020000  .0047-26.10115-43.09-22.22  .30 

no  40-100  .00  11  8.05  .67  .00  15  5  30 


69 


96030200  .0047-26.10115-43.09-22.39  .30  .0 

75  50-130  .00  9  8.05  .81  .00  20  10  25 


70 


960402:0fl  .00  47-26.10  115-43.09  -23.10 

160  65-170  .00  11  8.05  .44  .00  535 


.30 


.0 


71 


96053001  .0047-26.10115-43.09-17.77  .30 

125  90  -80  .09  11  8.05  .60  .00  3  0  25 
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96092604  .0047-26.10115-46.09-23.16  .30 

280  80  130  .10  10  8.05  .63  .00  33  15  30 


960926011 
55  15  SO 


.0 


.00  47-26.10  115-45.09  -23.18  .30 

.10  10  a. 05  .52  .00  IS  3  0  * 
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9S110101  .00  47-26.10  115-43.09  -26.00 

195  65  170  .DO  6  8.05  .42  .00  8  5180 


.30 


.0 


75 


96110102  .00  47-26.10  115-43.09  -26.00 

70  45-130  .00  10  8.05  .78  .00  23  10  25 


IJFDWN 

00  DILRTflTION 
+  +  COMPRESSION 


96110102 
IQQ  UQ  -SO 


UFDWN 

Oo  DILRTRTION 
COMPRESSION 


. 00  47-26 .10  11 5-y3 . 09  -26 . 00  .  SO 

DO  ID  a. 05  .ea  .do  is  s  lociK 
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96122001 
175  30  -90 


.00  47-26.10  115-43.09  -23.00  .30 

00  9  8.05  .87  .00  30  23  25 


+  +  COMPRESSION 


78 


96122001  .00  117-26.10  115-y5.09  -2S.00 

0  S5  -eo  .00  a  s.os  .79  .qo  20  is  iocsk 


t+  COMPRESSION 


STOZliiOO  .00  47-26.10  115-43.09  -24.03 
230  40  -60  .00  10  8.05  .75  .00  5  13  15 


UFDWN 

00  DILflTflTIGN 
+  +  COMPRESSION 


970311400  .00  47-26.10  115-43.09  -26.00 

235  50-100  .00  10  8.05  .85  .00  38  10  30 


UPDHN 

00  DILfITflTIGN 
+  +  COMPRESSION 


97052500 
240  40  I 


97052500  .00  47-26.10 

70  90-180  .10  10  a. 05 


97052500  .00  I17--26.10  115-y3.09  -24.39  .30  .0 

220  UQ  90  .20  10  9.05  .75  .00  10  S  5  « 


Appendix  B 


Fig.  B-1.  First  vibrational  mode  for  split-set  rock  bolt,  /  =  13  Hz. 


Fig.  B-2.  Third  vibrational  mode  for  split-set  rock  bolt,  fi  =  5\  Hz. 


Fig.  B-3.  Fifth  vibrational  mode  for  split-set  rock  bolt,  =  83.6  Hz. 


Fig.  B-4.  Sixth  vibrational  mode  for  split-set  rock  bolt,  7^  =  146.5  Hz. 
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Fig.  B-5.  Seventh  vibrational  mode  for  split-set  rock  bolt,  f,  =  232  Hz. 


Appendix  C 


Fig.  C-1.  10  kN  tension, /  =  28.6 /fz. 
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